Improved algorithms for splitting full matrix algebras 
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Abstract 

Let K be an algebraic number field of degree d and discriminant A over Q. Let A 
be an associative algebra over IK given by structure constants such that A = M n (K) 
holds for some positive integer n. Suppose that d, n and |A| are bounded. In a previous 
paper a polynomial time fF-algorithm was given to construct explicitly an isomorphism 
A — >• M n (K). 

Here we simplify and improve this algorithm in the cases n < 43, K = Q, and n = 2, 
with K = Q(V — 1) or IK = Q(\/— 3). The improvements are based on work by Y. Kitaoka 
and R. Coulangeon on tensor products of lattices. 
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1 Introduction 



x 



The following explicit isomorphism problem is important in computational representation the- 
ory: let IK be an algebraic number field, A an associative algebra over K. Suppose that A is 
isomorphic to the full matrix algebra M n (K). Construct explicitly an isomorphism A — > M n (K). 
Or, equivalently, give an irreducible A module. 

The algebra A is considered to be given by a collection of structure constants ^ijk G K. 
They form the multiplication table of A with respect to some IK basis a±, . . . , a m : the products 
didj can be expressed as 



7ijl«l + 7ij20 2 H h 7ijrrA 



In [5] a polynomial time ff-algorithm was given for the case of the problem, when n and 
the degree and the discriminant of K are all bounded. Applications were also outlined there, 
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including some parametrization problems of algebraic geometry. The methods of [5] are based 
on theoretical and algorithmic results on lattices and eventually boil down to enumerating short 
vectors in some lattices in real Euclidean spaces. 

In this paper we present considerable improvements of the methods of [5] in the case when 
the ground field K is the rationals, and n < 43. This is based on results of Kitaoka [6], [7] 
on tensor products of lattices, in particular we make use of the powerful result of Corollary 
[2J In Theorem [3] we prove a quantitative version, which allows further reduction in computing 
time. Some of Kitaoka's results have been extended by Coulangeon [3] from Q to imaginary 
quadratic fields. Using these we also obtain an improvement of the original algorithm when 
K = Q(v— 1) (Gaussian rationals) or K = Q(\/— 3) (Eisenstein numbers) for the case n = 2. 
The new algorithms are simpler and faster than the original ones. 

For the basic definitions and facts on lattices in real Euclidean spaces we refer to [2J, [9], 
and HQ]. 

2 Full matrix algebras over Q 

A (full) lattice L C M n is the free Abelian group generated by n linearly independent vectors 
bi, b2, . . . , b n . Then M = (t»ij t>2| • • • |b n ) is a matrix of L, and | det M\ is called the determinant 
of L, and is denoted by det L. We denote by X\(L) the Euclidean length of the shortest nonzero 
vector from L. The nth Hermite 's constant is 



where L is a full lattice in R n . Hermite proved that j n actually exists. The exact value of 7 n 
is known only for n G {1, 2, . . . , 8, 24}. 

We briefly recall now the definition of the tensor product of lattices, for more information 
see section 1.10 of Martinet [9], and section 7 in Kitaoka [TJ. Let L and M be two lattices in W 71 
and M n , respectively. The tensor product of Z modules L (gz M embeds in the straightforward 
way into M m ®m R n . This allows one to define L®Mas the set of integral linear combinations 
of the tensors x ® y from W 71 ®r MJ 1 where x G L and y G M. 

Note that, in terms of coordinates, L (g M can be viewed as the set (actually lattice) of 
m by n matrices over R which are integral linear combinations of dyads of the form xy T , 
where x G L and y G M. Note also that M m £*D M n is an Euclidean space with the law 
(xi eg) yi, X2 (g) y2) = (xi, X2)(yi, y2). In this setting the norm on the tensor product M™ 1 (g) M n 
is actually the Frobenius norm on the space of matrices M m ^ n (M.) . 

Let L be a full lattice in M. m . The dual L* of L consists of those vectors y G M m for which 
we have (x, y) G Z holds for every x G L. The supremum 7^ of Ai(L)Ai(L*) among full lattices 
L C lR n is the B erg e- Martinet constant, see PQ, [9]. It is known that 7^ < 7„ for every n. 

Following Y. Kitaoka [6] we say that a lattice L is of E-type if every minimal nonzero vector 
of L <g M is of the form x <g> y (x G L, y G M) for any lattice M. 

Kitaoka proved in [6], and in Theorem 7.1.1. of [7] the following. 

Theorem 1 (Kitaoka). If L is a lattice of rank at most 43, then L is of E-type. 

We remark here that by a theorem of Steinberg (see [TU], Chapter II, §9) for any integer 
n > 292 there exists a lattice L such that Xi(L (go L*) < Ai(L)Ai(L*). Thus, the conclusion of 
Kitaoka's theorem does not hold for larger values of n. 
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We can apply Kitaoka's theorem to maximal orders of the form A = QM n (Z)Q 1 C M n 
where Q G GL n 



Corollary 2. For any dimension n < 43 and subring A C M n (M) of the above form the nonzero 
matrices with the smallest Frobenius norm in A are of rank one. 

Proof. The order A is given by the transformation matrix Q. As in the proof of Theorem 1 
of [5], we see that as a lattice A = Q7J 1 <g> (QZ n )*. But Q17 1 is a rank n lattice, hence it is 
of -E-type, giving that the matrices of minimal norm in A are dyadic products of vectors from 
QUI 1 and (QZ* 71 )*, and hence have rank 1 as matrices from M n (R). □ 

Thus, when n < 43, then the smallest zero divisors in any A have rank one. In particular, 
no matrix from A of rank at least 2 can have minimal length. By using a modified variant 
of Kitaoka's original argument, we prove a slightly stronger, quantitative version of the latter 
statement. 

Recall that the rank of a nonzero tensor v G L <g> M is the smallest positive integer r such 
that v can be written as 

r 



v = ^x i <g>y i (1) 



i=l 

for some xi, . . . , x r G L and yi, . . . , y r G M. 

Theorem 3. Let L and M be lattices. Then for every tensor v G L <g> M of rank r we have 

|| v || > J^\ 1 (L®M). 

V % 

We remark that the above bound is meaningful when the dimension of L or M is at most 
43. Then obviously r < 43, and by Lemma 7.1.2 from [JJ for 2 < r < 43 we have 1 < r/7 2 . 
For r large the bound becomes trivial, because tends to zero as r grows. We shall need 
Lemma 7.1.3 from [7]: 

Lemma 4. Let A,B€L M n (R) be positive definite real symmetric matrices. Then we have 
Tr(AB) > n</detA^detB. 

Proof of Theorem O Let v G L <g> M be a tensor of rank r. Then v can be written in the 
form ([T]). Let L\ be the lattice generated by {xi, . . . ,x r }, and similarly let M\ be the lattice 
spanned by {yi, . . . , y r }. By the minimality of representation (JT]) the rank of these sublattices 
is r. Noting that 



Ml 2 



r 
i=l 



2 



^(x^Xj^y^yj) = Tr([(xi,x.,-)]^- =1 • [(yi,yj)]l j=1 ) 



and by using Lemma H] we obtain 

||v|| 2 >r(det[(x t ,x,)]-det[(y 4 , yj )]) 1/r . (2) 
Now let us assume for contradiction that ||v|| 2 < (r/7 2 )Ai(L (g) M) 2 . It follows that 

||v|| 2 < -tt(Ai(L)Ai(M)) 2 < -^(A 1 ( J Li)Ai(M 1 )) 2 . 
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Combining this with (j2J) we obtain 



7? (det[(x i ,x J )])^ (det^y,)])^-^ 7 - 

since [(xj, Xj)][ J=1 and [(y«, yj)][ J= i are Gram matrices for L\ and Mi, respectively. The con- 
tradiction finishes the proof. □ 

Using the known values of 7 r simple calculation gives that the minimal value of r/7! for 
2 < r < 8 is |, which is attained at r = 2. We remark that the bound of Theorem [3] is sharp, 
at least for r = 2. This is demonstrated by the hexagonal lattice A 2 < ffi 2 which is generated 
by the vectors 

Y 




It is known that A 2 is attains the Hermite constant, moreover this holds also for the dual lattice 
A 2 which is spanned by the vectors 




Some calculation shows that the minimal norm among the rank two tensors in A 2 (8> A* 2 is 
V2=^l-j= 3 = ^IX l (A 2 )X 1 (A*). 

In [5] it was shown that for the shortest nonzero matrix v G A = QM n (7j)Q~ l we have 
|| v ll < In, where 7™ is the Hermite constant. This can be strengthened as follows. Using again 
that A = QZ n ® (QZ n )*, we obtain that 

||v|| < Ai(Qz-)Ai((gz re )*) < y B , (3) 

where 7' is the Berge-Martinet constant. 



3 The modified IRS algorithm over Q for n < 43 

The input of the algorithm is an associative algebra A over Q given by structure constants. It 
is known that A is isomorphic to the full matrix algebra M n (Q). The objective of the algorithm 
is to find an element C G A which has rank one, when viewed as a matrix from M n (Q). 

The first four steps of the algorithm below are identical to the first four steps of the cor- 
responding algorithm from [5]. The last two steps of that method are replaced here by a new 
step 5: 

1. Construct a maximal order A in A. 

2. Compute an embedding of A into M n (R). This way we have a Frobenius norm on A. 
For X G A we can set ||X|| = ^/Tr(X T X). Also, via this embedding A can be viewed 
as a full lattice in M m , where m = n 2 . The length ||v|| of a lattice vector v is just the 
Frobenius norm of v as a matrix. 

3. Compute a rational approximation A of our basis B of A with a suitable precision. 
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4. Obtain a reduced basis b 1; . . . , b m of the lattice A C M. m by computing an LLL-reduced 
basis from A. The value c m of reducedness is ij m )~ (|) m 2 2~ . 

5. Generate all integral linear combinations 

C = ^ a ^ 

i=l 

where cti are integers, \a\\ < c m , until a C is found with rankC = 1. Output this C. 

Theorem 5. This algorithm is correct when n < 43. Moreover, it runs in ff-polynomial time. 

Proof. For the details and timing of steps 1-4 we refer to the proof of Theorem 1 from [5] . As 
a result of these computations we obtain a basis bi, . . . , b m of the lattice A such that 

||bi|| • ||b 2 || ■ ■ • ||b m || < c m ■ det(A) 

holds with 



m(rn — 1) 
2 



We recall the following bound by H. W. Lenstra [8]. 

Lemma 6. Let Y be a full lattice in IR m . Suppose that we have a basis b 1; . . . , b m ofT over Z 
such that 

||bi|| ■ ||b 2 || ■ •• ||b m || <c-det(T) 
holds for a real number c > 0. Suppose that 



v = ctjbj £ r, di £ Z. 



i=l 



T/ien we /iave lar» 1 < CtHt for i = 1, . . . , m. 

II II 

Let v £ A be a nonzero vector with minimal length. Then, on one hand, by Corollary [2] v is 
a matrix of rank one. On the other hand, when v is expressed as an integer linear combination 
v = Y^T=i a J°i then by Lemma i we have 

II v II 

I Oil < Cm j|bTj| - Cm - ( 4 ) 

This implies that there is indeed a rank one matrix C £ A among the linear combinations 
enumerated. 

To obtain the timing bound, we observe that at step 5 we enumerate at most (2c m + l) m 
linear combinations, and this value is bounded by our assumption n < 43. □ 

We remark that the upper bound \ai\ < c m which defines the domain to be searched at 
step 5 can be reduced. During the run of step 5 one may update the quantity d which is the 
actual minimum of the numbers (7;?/^) Il^ll over the matrices C enumerated so far (here r is 
the rank of C). From 01]), Theorem [31 and ([3]) it follows that 

M < c m ^^-. (5) 
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4 Two by two matrices over imaginary quadratic fields 

Here we consider possible extensions of the improvements obtained over Q to other number 
fields. In general not much is known about tensor products of lattices over general number 
fields. On the positive side, Coulangeon j3] extended some of Kitaoka's results to imaginary 
quadratic fields. 

Let K denote an imaginary quadratic number field Q(V — d) where d is a square-free positive 
integer. By O we denote the ring of algebraic integers in K: O = Zl + l,y / —d + Z 1+ ^~^ if 
d = — 1 modulo 4 and O = Zl + Zyf—d otherwise. In the next discussion we consider K (and 
hence O) to be embedded into C. 

Let A be a central simple algebra of dimension 4 over K isomorphic to M 2 (K) and let A be 
a maximal order in A. We assume that we are given an embedding <fi of A into M 2 (C). From 
the theory of central simple algebras over number fields we know (see Corollary 27.6 in Reiner 
[TT] ) that there exists a matrix B G M 2 (C) and a fractional ideal / of O such that 

In other words, there exists a full CMattice L in C 2 (a finitely generated (9-submodule of C 2 
that spans C 2 as a linear space over C) such that 

A = {A G M 2 (C) : AL C L} . 

In fact, for our specific A, we can take L = B 

Let (, ) stand for the standard Hermitian bilinear form on C 2 . The linear extension of the 
mapping u <g) v \-> A UjV where A UjV w = (v, w)u gives an identification of M 2 (C) with C 2 (g>c C 2 . 
By this identification, the tensor square of the standard Euclidean norm of C 2 becomes the 
Frobenius norm of matrices. Also, A is identified with L ®o L* where 

L* = [u G C 2 : (u, v) G L for every 

and L ®o L* is just the additive subgroup of C 2 ®c C 2 spanned by the tensors of the form u <g> v 
where u G L and 

Remark. We can speak about the rank of an element x G L ®o L* in two ways. One is the 
minimal positive integer r such that x can be written as 

r 

x = ^Xi<g>yi, (6) 
i=i 

for some vectors Xj G L and y« G L*. The other possible notion of rank is the rank of x as 
a matrix from M 2 (C). The two notions are not the same^. As an exampl^], let d — 5, and 
A = M 2 (0) and consider the matrix 




: It is not hard to show that the two notions of rank coincide if O is a principal ideal ring. 
2 We thank Geza Kos for suggesting this example. 
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from A. We have det C = 0, hence C has rank 1 as a matrix from M 2 (C). Moreover, using the 
fact that 3 and 2 are irreducible elements in O, we see that C is not a decomposable tensor 
from O 2 ® O 2 . 

We shall use the term rank in the former sense. We note also, that in the minimal repre- 
sentation ([6]) the vectors Xj and y, are linearly independent over K. For a proof we refer to 
Lemma 3.1 in [3]. 

Let M be an (9-submodule of C 2 generated by two linearly independent vectors. The 
determinant of M is defined as the following Gram matrix 

detM = det(< Vl ' Vl > j Vl ' V2 l 

V( V 2, v l) (V 2 ,V 2 ) 

where vi, v 2 is any basis for M. 

Following the notation of [3], we denote by jh (M) the quantity 

||v||7(det M)i, 

where v is a shortest nonzero vector from M. 

Let D be the discriminant of IK (we have D — d if d = 3 modulo 4 and D = 4d otherwise.) 
It is known (last paragraph of Subsection 2.1 in [3]) that 7ft (M) = r y(M)vD/2, and hence 



lh (M) = 7(M)v^/2 < 74V^/2 = y/D/2, (7) 

where 7(M) is the ratio of ||v|| 2 and the fourth root of the determinant of M, considered as a 
Z-lattice of rank 4, and 74 = v2 is the Hermite constant. 

Let us define r(A) as the ratio between the squared length of the shortest rank 1 element 
of 0(A) and that of the shortest rank 2 element in 0(A). 

Proposition 7. We have 

r(A) < ^(Mhft(M'), 

where M is an O-sublattice of L generated by two linearly independent vectors over IK and and 
M' is an O-sublattice of L* generated by two linearly independent vectors over K. 

Proof. Indeed, let v and w be shortest nonzero vectors from L and L*, respectively. Also, let 
w be a shortest nonzero vector of rank 2 from A = L ®o L*. We have then 



r(A) = 

Similarly to the rational case ()2]) we obtain 



|v|| 2 • |w|| 2 

IMI 2 ' 



IMI 2 > 2(detM) 1/2 (detM') 1/2 

for some sublattices M < L and M' < L* which are spanned by two linearly independent 
vectors over IK (see Proposition 3.2 from [3] for the details). 

Let v' and w' be shortest nonzero vectors from M and M', respectively. Clearly we have 
|| v || < ll v 'l|; an d ll w ll ^ ll w '||- By putting all these together we obtain 

||v|| 2 -|w|| 2 llv'11 2 ■ Iw'll 2 1 

r(A) = HF - 2(detM)V2( de tM0V2 = 2^^)- 

□ 
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For d = 1 by ([7]) and Proposition [7] we have r(A) < |\/2 ■ \/2 = 1. Similarly, for d = 3 we 



find that 



We have obtained the following: 



/A , 1/3/3 3 
r ' A » £ 2V2V2 = 4 <L 



Proposition 8. For d = 1, at least one of the smallest element of 0(A) with respect to the 
Frobenius norm has rank one. For d = 3 every smallest element of 0(A) has rank one. □ 

The following example shows that over the Gaussian rationals it does indeed occur that a 
shortest nonzero element of A has rank 2. Let d = 1 let L be (9-sub module of C 2 generated by 
(1, 0) T and (^=, ^) T and let A = {A G M 2 (C) : AL C L}. Then 



A 



1 + i/a 6 

c e 



a, b,c,e G 0, 
a + c = a + 6 = 6 + e = c + e = mod (1 + z) 
a + 6 + c + e = mod 2 



and the identity matrix is one of the elements of A having the smallest Frobenius norm. 

Next we outline a direct, elementary proof of inequality (JTJ. Let K = Q(\/— d) be an 
imaginary quadratic number field, O be the maximal order of K. Suppose further that O is a 
principal ideal ring. 

Lemma 9. Let z be any complex number. Then there exists an element a G O such that 
\z — al < K, where 



K 



g2j d = 3 mod 4, 



/a+T 



otherwise. 



The proof is a simple argument from elementary geometry which we omit here. For d = 1, 



we have k = for d = 2, k 



2 ' 



2 ' 



for d = 3, k 



3 ' 



for d = 7, ft = ^y^, for d 



11, 



k = -4^. In these cases k < 1. For for d = 5, 6, 10 and for d > 11 we have k > 1. Let us define 
r = + lj. Then obviously we have k/t < 1. 

Proposition 10. Suppose that k < 1 holds, and let M be an O-submodule of C 2 generated by 
two linearly independent vectors over K. T/ien 



7/x(M) < 



r 



.2 



Proof. Let v, w be a basis of M such that v is a shortest nonzero vector from M. Such a basis 
exists because O is a principal ideal ring. Apply now the preceding lemma for z = r-^^y. 
There exists an a G O be such that 



(v, w) 



(v,v) 



< 



K 



T 



and put w' = w — -v. Then rw' G M and hence llw'11 > — . Furthermore 

J- 7- II — T 

|(v,w')| = 



a 

-v) 

r 



(v, w) (v, v) < -(v, v). 

r t 



S 



The Gram determinant does not change if we switch from the basis v, w to v, w', hence 
detM = det (>; V v > } { ( ;' W j ) )=(v,v)(wy)-|(v,w')P 

> |i-£lllvl" 



□ 



For d — 1 the Proposition gives 'yh(M) < \[2 = WD~/2. For d = 2 we obtain 7h(M) < 
2 = y/D/2. For d = 3 our bound is 7h(M) < ^ = y/D/2. For d = 7 the proposition gives 

7/j(M) < a/| < \/| = yf ■ ^or these values of d the ring (9 is a principal ideal ring, hence we 
have proved (j7|). 

The improved algorithm when c£ = 1 or d = 3 

We can achieve an improvement of the algorithm of Section 3 from [5] for n = 2 and d — 1 
or g? = 3, i.e. for the case of two by two matrices over the Gaussian rationals or over the 
Eisenstein rationals. These cases of the explicit isomorphism problem occur when one considers 
parametrization of Del Pezzo surfaces of degree 8, see Section 4 in [3]. Our method may present 
a viable alternative there to solving norm equations. 

Our improvement over the method of [5] is very similar to that of the algorithm over Q. 
Suppose therefore that K is either Q(a/— 1) or Q(v— 3), and we have as input an algebra A 
over IK specified by structure constants. We assume that A = M 2 (K). The first four steps 
of the method in [5] construct a maximal order A, an embedding (j) : A — > M2(C), a Z linear 
embedding $ of A into 1R 8 which maps an y G A to 

$(y):=(3?0(y),S0(y))GM 8 . 

The image $(A) is a full Z lattice in M 8 . Note that the (real) Euclidean norm ||$(?/)|| is the 
same as the Frobenius norm ||0(y)|| inherited from M 2 (C). 

Moreover, the first four steps return a Z basis bi, . . . , b§ of $(A) for which we have 

||bi|| ■ |b 2 || ■ ■ ■ ||b 8 || < c 8 det$(A). 

From thispoint on we can simplify the search for a zero divisor. By Proposition [8] it suffices 
to enumerate^ the elements of $(A) which have minimal norm until a zero divisor in A is found. 
To this end, we generate all integral linear combinations 

8 

v = y^a^ba, 

i=l 

where ctj are integers, |a»| < eg, until a v is found for which the matrix y G A with $(?/) = v 
is of rank one. The bound on the integers «j follows from Lemma El like in the rational case. 
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3 For d = 3 it suffices to find just one element with minimal norm. 
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